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Abstract. We present an astrometric analysis of the binary systems AB Dor A / AB Dor C and 
AB Dor Ba /AB Dor Bb. These two systems of well-known late-type stars are gravitationally associated and they 
constitute the quadruple AB Doradus system. From the astrometric data available at different wavelengths, we 
report: (i) a determination of the orbit of ABDorC, the very low mass companion to ABDorA, which confirms 
the mass estimate of Q.O9OM0 reported in previous works; (ii) a measurement of the parallax of ABDorBa, which 
unambiguously confirms the long-suspected physical association between this star and ABDorA; and (iii) evi- 
dence of orbital motion of ABDorBa around ABDorA, which places an upper bound of 0.4 M© on the mass of 
the pair AB Dor Ba / AB Dor Bb (50% probability). Further astrometric monitoring of the system at all possible 
wavelengths would determine with extraordinary precision the dynamical mass of its four components. 

Key words, astrometry - stars: kinematics - stars: binary - stars: low mass stars: individual (AB Doradus) - stars: 
individual (Rst 137 B) - 



1. Introduction 

■ Astrometry, one of the most classical astronomical disci- 
, plines, provides unambiguous mass estimates of celestial 
' bodies via observations of the orbits of binary or multi- 
, pie systems, as projected on the sky (see, e.g., Kovalevsky 
1995). The precise determination of stellar masses is fun- 
damental in astronomy, as this parameter is the primary 
input to test stellar evolutionary models that provide 
widely used mass-luminosity relations. In particular, the 
calibration of the mass-luminosity relation for the lower 
end of the main sequence is of special interest, since it 
permits the derivation of the physical properties of very- 
low-mass (VLM) stars and substellar objects. However, a 
model-independent measurement of the mass of these ob- 
jects is a most demanding task that requires the astromet- 
ric follow-up of VLM stars in binary systems (e.g. Lane et 
al. 2001; Bouy et al. 2004; Golimowski et al. 2004; Close 
et al. 2005). One of the few VLM objects with astrometric 
detections is ABDorC, the companion to ABDorA. 
ABDorA (=HD 36705) is an active Kl star only 14.9 pc 
away from the Sun. Due to its ultrafast rotation (0.514 
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days; Innis et al. 1985), ABDorA is a strong emitter 
at all wavelengths, and it has been extensively observed 
from radio to X-rays (Lim et al. 1992; Mewe et al. 1996; 
Vilhu et al. 1998; Giidel et al. 2001). ABDorA possesses 
a low-mass companion, ABDorC, which induces a reflex 
motion first detected by very-long-baseline-interferometry 
(VLBI) and the Hipparcos satellite (Guirado et al. 1997). 
Recently, Close et al. (2005) [CLG] obtained a near- 
infrared image of ABDorC, providing the first dynam- 
ical calibration of the mass-luminosity relation for low 
mass, young objects. AB Dor A has another physical com- 
panion, ABDorB (=Rossiter 137B, =Rstl37B), a dM4e 
star, which is also a rapid rotator with a 0.38 day pe- 
riod and is separated from ABDorA by 9" (Lim 1993). 
Based on their young age (CLG), common proper mo- 
tions, and common radial velocities (Innis, Thompson & 
Coates 1986), it is believed that both stars may be associ- 
ated. In turn, CLG found ABDorB to be a tight binary 
(ABDorB^ABDorBa and ABDorBb). 
AB Dor C is the first calibration point for evolutionary 
tracks in the young VLM regime. From comparison with 
theoretical predictions, CLG found that the dynamical 
mass of AB Dor C is almost twice than predicted by evolu- 
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tionary models (Chabrier et al. 2000), which suggests that 
models tend to underpredict the mass of young VLM ob- 
jects. In this context, a precise estimate of the dynamical 
mass of ABDorC is extremely important. In this paper 
we report the details of an improved method to deter- 
mine the mass of ABDorC, which confirms the value of 
0.090 M0 given by CLG. We also report on the sky motion 
of AB Dor Ba, which shows a nearly-identical parallax to 
that of AB Dor A and evidence of the long-term orbital 
motion around AB Dor A. 



2. Astrometric Data 

In Table 1 we summarize the available astrometric 
data of the AB Doradus system, which include abso- 
lute positions of ABDorA and ABDorBa, relative po- 
sitions of the 9" pair AB Dor A /AB Dor Bb, and rela- 
tive positions of the closer pairs AB Dor A / AB Dor C 
and ABDorBa/ABDorBb. New absolute positions of 
AB Dor Ba are presented in this table; they have been ob- 
tained from the same VLBI observations that were used 
to make the astrometric analysis of AB Dor A reported by 
Guirado et al. (1997). Given the 9" separation, ABDorA 
and AB Dor Ba lie within the primary beam of each of 
the telescopes and thus can be observed simultaneously 
for efficient cancellation of atmospheric systematic errors. 
The interferometric array has much finer resolution (a few 
milliarcseconds) and, therefore, the interferometric data 
for AB Dor Ba could be extracted and processed following 
the same procedures as described in Sect. 2 of Guirado et 
al. (1997) for ABDorA. This in-beam technique is widely 
used in VLBI observations (e.g. Marcaide & Shapiro 1983; 
Fomalont et al. 1999). On the other hand, the relatively 
low brightness of ABDorBa {V=12.6; Collier Cameron 
& Foing 1997) explains the absence of Hipparcos data for 
this star. In Sect. 3, we revisit the astrometry of the dif- 
ferent pairs shown in Table 1. 



3. Astrometric Analysis 

3.1. ABDorA /ABDorC: Orbit Determination 

The infrared image of AB Dor C provided the astromet- 
ric data that was used by CLG to constrain the elements 
of the reflex orbit. The weakness of this procedure was 
that the relative position AB Dor A / AB Dor C was not in- 
cluded in the flt, rather it was only used as a discriminator 
of the orbits that plausibly fit the VLBI/Hipparcos data. 
In this section, we re-cstimatc the mass of AB Dor C using 
a much improved method that estimates the reflex orbit 
of AB Dor A by simultaneously combining both the exist- 
ing VLBI/Hipparcos ABDorA astrometric data and the 
near- infrared relative position of AB Dor C. Following the 
classical approach, we modeled the (absolute) position of 
AB Dor A (a, 6) at epoch t from the expressions: 



+ Sa{t, Xi, X2, X3, X4, P, e. To) 
5{t)=S{to) + iis{t-to)+T:Ps 

+ Ss{t, Xi, X2, Xs, X4, P,e,To) 



(1) 



where to is the reference epoch, Ha, Us are the proper 
motions in each coordinate, tt is the parallax, P^ and Pg 
are the parallax factors (e.g. Green 1985), and Sa and 
Ss are the reflex orbital motions in a and 5, respectively. 
The astrometric parameters {a{t{)), 5(to), ^^a, ^^S, and 
tt) are linear in Eq. (1). The reflex motion, Sa and Ss, 
depends on the seven orbital parameters, namely, a, 
i, lJ, f7, P, e, and Tq. Wc have used the Tliieles-Imies 
coefficients (Green 1985), represented by Xi, X2, X3, X4, 
which are defined as combinations of a, i, u), Q.. These 
coefficients behave linearly in Eq. (1), leaving only three 
non-linear parameters (P, e, and Tq) to solve for in our 
weighted-least-squares approach. 

Since our fitting procedure estimates the orbital parame- 
ters of the reflex motion of ABDorA, the relative sepa- 
ration AB Dor A / AB Dor C provided by the infrared data 
(Aq', A5') at epoch t' is included in the fit via the cor- 
responding orbital position of the primary star according 
to the definition of the center of mass of the system: 
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where q is the mass ratio mc/nia, with rria being the mass 
of the primary and rric the mass of the companion. The 
combination of data types in the same fit is reflected in 
the definition of the to be minimized: 
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a{t) = a{to) + Hait- to) + Pa 



where the tr's are the corresponding standard deviations 
(Table 1) and the circumflexed quantities are the theoret- 
ical values of our a priori model. The virtue of the defi- 
nition of in Eq. (3) is that the linearity of the orbital 
parameters is conserved as long as the mass ratio q is not 
adjusted in the fit. In consequence, rria is a fixed param- 
eter in our fit (we adopted the value of O.865±O.O34M0, 
as given by CLG). The mass of the secondary (mc) will 
be estimated via an iterative procedure that we outline 
below: 

1. We set a priori values of the three non-linear param- 
eters (P, e, and To). In particular, we sample the fol- 
lowing parameter space: < P < 30 years, 
1990.0 < To < 1990.0 -hP, and 0< e < 1. To start 
this procedure, we need an initial value of rric. We take 
a value of 0.095 M©, which corresponds to the central 
value of the rric interval given in Guirado et al. (1997). 
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AB Dor A 



Epoch 


Instrument 


q;(J2000) 


<5(J2000) 


Reference 


1990.3888 


Hipparcos 


S** 28™ 44f 77474 ± 0f00026 


-65° 26' 56"2416 ± 0"0007 


(1) 


1990.5640 


Hipparcos 


5''28'"44f 78652 ± Of 00025 


-65° 26' 56''2272 ± 0''0007 


(1) 


1991.0490 


Hipparcos 


28"" 44f 77578 ± Of 00024 


-65° 26' 56"2615 ± 0"0007 


(1) 


1991.5330 


Hipparcos 


g;, 28™ 44f 78942 ± Of 00025 


-65° 26' 56''0757 ± 0''0008 


(1) 


1992.0180 


Hipparcos 


28™ 44f 78202 ± Of 00024 


-65° 26'56"1160 ± 0"0009 


(1) 


1992.2329 


VLBI 


28™ 44f 77687 ± Of 00019 


-65° 26' 56"0049 ± 0"0007 


(1) 


1992.6849 


VLBI 


5'" 28™ 44f 80124 ± Of 00018 


-65° 26' 55''9395 ± 0"0006 


(1) 


1993.1233 


VLBI 


gh. 28™ 44f 78492 ± Of 00024 


-65° 26' 55"9137 ± 0"0008 


(1) 


1994.8137 


VLBI 


28™ 44f 81768 ± Of 00019 


-65° 26' 55"6866 ± 0"0005 


(1) 


1995.1425 


VLBI 


2gm 80247 ± Of 00027 


-65° 26'55"6248 ± 0"0011 


(1) 


1996.1507 


VLBI 


gh 28"" 44f 81137 ± Of 00013 


-65° 26' 55"4852 ± 0"0003 


(1) 


1996.3607 


VLBI 


5''28'"44f81776 ± Of 00018 


-65° 26' 55:'3785 ± 0:'0010 


(1) 






ABDorBa (=Rstl37B) 




Epoch 


Instrument 


a(J2000) 


5(J2000) 


Reference 


1992.2329 


VLBI 


gft 28™ 44f 39520 ± Of 0007 


-65° 26' 47"0676 ± 0"0024 


(2) 


1992.6849 


VLBI 


gh 28™ 44f 41973 ± Of 0006 


-65° 26' 47" 0047 ± 0"0021 


(2) 


1993.1233 


VLBI 


gh 28™ 44f 40441 ± Of 0008 


-65° 26' 46"9869 ± 0"0028 


(2) 


1994.8137 


VLBI 


-h 2grn 44f43(587 ± of 0007 


-65° 26'46"5528 ± 0"0018 


(2) 


199().ir)()7 


\'LB1 


."•)'' 28"' ii: 12812 ± ().^()()()5 


-()5-' 2()' 1():':)77:! ^ ()''0010 


(2) 






Relative Position AB Dor A 


- AB Dor Ba 




Epoch 


Instrument 


Separation 


RA. (°) 


Reference 


1929 


_ 


10"0 


339 


(3) 


1985.7 


AAT 


9:'3 ± 0''3 


344 ± 5 


(4) 


1993.84 


ATCA 


8"90 ± 0"02 


345.2 ± 0.1 


(5) 


1994.2 


Dutch/ESO 


8!'9 ± 0''1 


344.7 ± 0.3 


(6) 


2004.093 


VLT/NACO 


9!'01 ± Q'.'Ol 


345.9 ± 0.3 


(7) 






Relative Position ABDorA 


- ABDorC 




Epoch 


Instrument 


Separation 


P.A. (°) 


Reference 


2004.093 


VLT/NACO 


0"156 ± 0"010 


127 ± 1° 


(7) 






Relative Position ABDorBa 


- AB Dor Bb 




Epoch 


Instrument 


Separation 


P.A. (°) 


Reference 


2004.098 


VLT/NACO 


0"062 ± 0"003 


236.4 ± 3.33° 


(8) 



(1) Guirado ct al. (1997); (2) this paper; (3) Jcffcrs ct al. (1963); (4) Innis et al. (1986); (5) J. Lim, 
personal communication (6) Martin & Brandner (1995); (7) Close et al. (2005); (8) Brandner et al. 
in preparation 



2. To find a minimum of x^i s-s defined in Eq. (3), we 
used an iterative method, based on the Brent algo- 
rithm (Press et al. 1992). The minimum is defined such 
that the difference between the rc}duc;c;d-x^ of two suc- 
cessive iterations is significantly less than unity. 

3. From the resulting orbital parameters, wo then use 
Kepler's third law [m^/(mo + rric)^ = af/P^, with ai 
the semimajor axis of the reflex orbit] to estimate the 
mass rjic of AB Dor C. 

4. We iterate the least squares fit (step 2) using as a priori 
values the new set of adjusted orbital parameters, and 
estimated rric. 

5. A final set of orbital parameters is obtained once the 
estimated nic is self-consistent, that is, the difference 
between the value of rric calculated in step 3 from con- 
secutive sets of adjusted orbital parameters is negligi- 
ble (i.e. <O.OO1M0). 

The resulting orbital parameters, and the estimate of the 
mass of AB Dor C, are shown in Table 2 and represented in 



Fig. 1 . These values are fully compatible with those given 
by CLG. However, our method shows the robustness of 
the determined orbit. Despite the wide range of param- 
eter space investigated, the solution found for the reflex 
orbit of AB Dor A is unique (see Fig. 2). This is a remark- 
able result: for astrometrically determined orbits, Black 
& Scargle (1982) predicted a coupling between the proper 
motion and the orbital wobble, resulting in an underesti- 
mation of the period and semimajor axis. This coupling 
is present in the VLBI/Hipparcos data that only covers 
51% of the reflex orbit. Our least-squares approach copes 
partially with this efli'ect, since the astrometric and or- 
bital parameters are estimated simultaneously. However, 
the VLT/NACO data not only extends signiflcantly the 
observing time baseline, but represents purely orbital in- 
formation, independent of proper motion and parallax ef- 
fects. In practice, the combination of astrometric data of 
the primary star with astrometric data of the relative orbit 
improves the fit dramatically, constraining the orbital pe- 
riods allowed by the astrometric data of the primary only. 
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Table 2. J2000.0 astrometric and orbital parameters of 
AB Dor A 




Fig. 1. Above: orbits of the pair ABDorA (inner el- 
lipse) and ABDorC (outer ellipse). Below: blow up of the 
dotted square in the figure above. VLBI and Hipparcos 
data points are marked in ABDorA's orbit, while the 
VLT/NACO ABDorC position relative to ABDorA is 
indieated in ABDorC's ellipse. The star symbols over 
the orbits correspond the astrometric predictions at epoch 
2004.093, based on the orbital elements given in Table 2. 



Parameter 









— fi^ 9fi 'I'l Q'i'i 


1 1 ( ^ \!"r 


n nn77 + n nnn? 


1 1 r ( 9\ TPtJf^P \TV ^ \ ' 
f-tiO I Oil y i. J. 




n (arcsec): 


0.0664 ± 0.0005 


P(yr): 


11.76 ±0.15 


ffli (arcsec): 


0.0322 ± 0.0002 


e: 


0.60 ± 0.04 


i(deg): 


67 ±4 


aj(deg): 


109 ±9 


O(deg): 


133 ±2 


To-. 


1991.90 ±0.04 


mc (M©)': 


0.090 ± 0.003 



" The reference epoch is 1993.0. Units of right ascension are 
hours, minutes, and seconds, and units of dechnation are de- 
grees, arcminutcs, and arcsoconds. 

Mass range obtained from the period and semimajor axis 
via Kepler's third law. The mass adopted for the central star 
ABDorA was 0.865±0.034 Mq. 



■■■■■■r-:^r-,t- 



In our case, the constraint is such that the only allowed 
period is 11.76±0.15yr. In general, our results show the 
combination of different techniques is more effective than 
any one technique alone. 

3.1.1. Error Analysis 

Our least-squares procedure provides formal errors for the 
adjusted orbital parameters. However, other systematic 
contributions need to be taken into account. In particu- 
lar, this includes the uncertainty associated to the mass 
of ABDorA, which is a fixed parameter in our analysis 



Mass (Mq) 

Fig. 2. Result of the exploration of the ABDorA reflex 
orbit. A well-defined minimum is found for a mass com- 
panion of 0.090 Mq. See Sect. 3.1. 



(O.865±O.O34M0). To estimate this error contribution, we 
altered the mass of AB Dor A by one standard deviation 
and repeated the fitting procedure to obtain the change in 
the orbital parameters and the mass of AB Dor C. We note 
that this a conservative approach, since this technique fails 
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Table 3. J2000.0 VLBI astrometric parameters of 
AB DorBa 



Parameter 


a 


5 28 44.4123 ± 0.0002 


S": 


-65 26 46.9974 ±0.0015 


IXa (syr"^): 


0.0085 ± 0.0002 


(arcsecyr""*^): 


0.134 ± 0.0012 


TT (arcscc): 


0.0666 ± 0.0015 



" The reference epoch is 1993.0. Units of right ascension are 
hours, minutes, and seconds, and units of dechnation are de- 
grees, arcminutes, and arcseconds. 

to account for the correlation of ma with the rest of the 
parameters. The resulting parameter changes were added 
in quadrature with the formal errors of our fit (see Table 
2). As expected, the 0.003 Mq standard deviation of the 
mass of AB Dor C is dominated by the uncertainty in the 
mass of ABDoiA, while the standard deviations of the 
rest of the parameters are dominated by the statistical er- 
rors. 

We also checked the dependence of our results on the 
choice of the a priori value of rric in step 1 of our fit- 
ting procedure (Sect. 3.1). We found that the results are 
insensitive to this choice. The postfit residuals of the po- 
sitions of ABDorA exhibits an rms of ~lmas at each 
coordinate, consistent with the standard errors, and with 
no evidence, within uncertainties, of any further orbiting 
companion to AB Dor A. 

3.2. VLBI Astrometric Parameters of AB DorBa 

Innis et al. (1985) presented radial velocity measurements 
of AB DorBa, the 9" companion to AB Dor A. Their mea- 
surements do not differ from those of AB Dor A within the 
uncertainties. Additionally, Innis et al (1986) and Martin 
& Brandner (1995) reported close agreement between the 
proper motions of both stars. These results are strong ar- 
guments in favor of a physical association of both stars. We 
used the VLBI (absolute) positions of AB Dor Ba given in 
Table 1 to derive the parallax and proper motion via a 
least-squares fit. The results of this fit are presented in 
Table 3, which shows that the parallax of AB DorBa is 
coincident with that of AB Dor A to within the uncertain- 
ties, which provides independent and conclusive evidence 
for the association of both stars. Comparison of Table 1 
and Table 3 shows that the proper motion of AB Dor Ba 
derived from the radio data appears significantly differ- 
ent to that of ABDorA. Given the relatively small un- 
certainty of our determination, this does not contradict 
previous (and coarser) measurements of common proper 
motion. Rather, we interpret this extra proper motion of 
AB DorBa towards the south-east as a result of the or- 
bital motion around ABDorA (see Sect. 3.3). 
The postfit residuals of AB Dor Ba show a systematic sig- 
nature, both in right ascension and declination, which 
corresponds to a relatively high rms of ~4mas. The 



short time span between our separate VLBI observations 
makes it unlikely that this signature is an effect of the 
long-term gravitational interaction of AB Dor Ba with 
ABDorA. Rather, this signature could be assigned to the 
0.070" companion (ABDorBb) of AB DorBa seen in the 
VLT/NACO observations reported by CLG. As for the re- 
vision of the reflex orbit of AB Dor A, we attempted to get 
estimates of the orbital elements of the reflex motion of 
AB Dor Ba by combining the radio data with the VLT rel- 
ative position between AB DorBa /ABDorBb (Table 1). 
However, our analysis did not yield useful bounds to the 
mass of this pair, showing that the number of data points 
is still insufficient and, more likely, they do not properly 
sample the expected short period of this tight pair. 

3.3. ABDorA /AB DorBa: evidence of orbital motion 
of AB DorBa 

As stated in the previous section, evidence of the motion 
of AB Dor Ba around AB Dor A can be obtained from the 
radio data alone. In order to get more precise evidence 
of this orbital motion, we augmented our data set with 
relative positions ABDor A/ABDorBa found in the 
literature (see Table 1). We then corrected all relative 
positions ABDor A/ ABDor Ba for the reflex orbital 
motion of ABDorA (Table 2), effectively referring the 
positions of AB Dor Ba to the center of mass of the 
ABDor A/ AB Dor C system. 

We attempted to constrain the relative orbit of 

AB Dor A / AB Dor Ba following a similar analysis to 
that described in Sect. 3.1, fltting only the 7 parameters 
of the relative orbit. We sampled all possible periods up 
to 5000 years and eccentricities from to 1. We selected 
as plausible orbits those whose reduced-x'^ differs by 25% 
of the minimum. For each plausible orbit, the mass of the 
complete system was estimated from Kepler's third law, 
now expressed in terms of the parameters of the relative 
orbit: 

(a/7r)^ ,^ ,^ 

= M(A+C) + M^Ba+Bb) (4) 

where Af(^_(.(7) and M(^Ba+Bb) are the combined masses of 
AB Dor A/AB Dor C and AB Dor Ba / AB Dor Bb, respec- 
tively, a is the relative semimajor axis (arcsec), tt is the 
parallax (arcsec; Table 2), and P is the period (yr). The 
poor coverage of the orbit favors a correlation between the 
orientation angles and the eccentricity, allowing a wide 
range of orbital parameters that flt our data equally well. 
However, a similar correlation between P and a imposes a 
constraint on the determination of the mass of the system 
via Eq. (4), which is represented in the histogram of Fig. 
3. From the plausible orbits selected, more than 50% cor- 
respond to a total mass of the AB Doradus system in the 
interval 0.95— 1.35 (see Fig. 4 for examples of plausible 
orbits). Larger masses are not excluded, but the required 
orbital configurations for masses outside this range occur 
with significantly reduced probability. 
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Mess of the ABDor system (Msun) 



Fig. 3. Histogram of plausible orbits for the relative orbit 
of ABDorBa around ABDor A. More than 50% of the 
plausible orbits correspond to a total mass of the system 
in the range 0.95-1.35 Mq. 



If wc assume the total mass of the AB Doradus sys- 
tem lies in the interval 0.95— 1.35Mq, the combination 
with our estimate of Afj^+c) (O.956±O.O35M0; see Sect. 
1) suggests an upper bound to the mass of the pair 
ABDorBa/ABDorBb of 0.4 Mq. This upper limit to 
M( Ba+ Bb) looks too coarse to calibrate evolutionary mod- 
els. Nevertheless, it can be transformed into a bound to the 
age of this pair. To do this, we used the JC-band 2MASS 
photometry of AB Dor Ba, and the /T-band difference be- 
tween ABDorBa and ABDorBb reported by CLG. The 
comparison with Baraffe et al. (1998) isochrones suggests 
an age for this pair in the range of 50— 120Myr. This range 
is compatible with previous values of the age of AB Dor Ba 
(30-100 Myr; Collier Cameron & Foing 1997). However, 
our age estimate for AB Dor Ba /AB Dor Bb is not conclu- 
sive: first, the masses of the individual components arc yet 
to be determined, and second, there are indications that 
the evolutionary models might need revision, since they 
tend to underprcdict masses for very young objects below 
0.3 Mq (CLG; Reiners et al. 2005). 





-2 -2.5 -3 -3.5 -4 

RA (") 

Fig. 4. Above: positions of ABDorBa with respect to 
the center of mass of ABDor A/AB Dor C (see Table 1) 
and several allowed orbital solutions. The displayed or- 
bits correspond to a total mass of the system in the 
range 0.95-1.35 Mq with periods of 1400, 2300, and 4300 
years. The cross at the origin indicates the position of 
ABDor A/ ABDor C. Below: blow up of the region con- 
taining the measurements. 



4. Summary 

We have revisited the different orbits in the quadruple 
system in ABDoradus. Paradoxically, this system, where 
the measurement of precise radial velocities is difficult 
due to the fast rotation of the main components, has be- 
come an extraordinary target for astrometric techniques 
in different bands of the electromagnetic spectrum. From 
our analysis of the available data, we have re-estimated 
the mass of the VLM star ABDorC by using a least- 
square approach that combines the data from radio, opti- 
cal, and infrared bands. Although the data do not cover 
a full orbit, the mass and orbital elements of ABDorC 
are strongly constrained and fully compatible with those 
reported by CLG. Further monitoring of the reflex orbit 



of ABDor A via VLBI observations, and of the relative 
orbit ABDor A /ABDorC via VLT/NACO observations, 
will result in independent estimates of the masses of the 
components of this pair. From the absolute radio posi- 
tions of AB Dor Ba, we have determined the absolute sky 
motion (i.e. not referred to the motion of ABDor A) of 
this star and, in particular, its parallax, which is identical, 
within the uncertainties, to that of AB Dor A. This con- 
firms the association of both stars. The mass of AB Dor C 
serves as a precise calibration point for mass-luminosity 
relations of young VLM stars. Likewise, other compo- 
nents of AB Doradus may provide new calibration points 
for slightly higher masses. We have found evidence for 
the long-term orbital motion of AB Dor Ba /AB Dor Bb 
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around ABDor A/ABDorC. From an exploration of the 
multiple orbits that fit the available data we find that 
the most probable mass upper limit of the pair is 0.4 Mq. 
This limit maps into an age range of 50—120 Myr using the 
isochrones provided by Baraffe et al. (1998). Further mon- 
itoring with the appropriate sampling, both in radio and 
infrared, should provide the orbital elements of both the 
relative and reflex orbits of the pairs AB Dor A / AB Dor C 
and AB Dor Ba / AB Dor Bb, from which would follow pre- 
cise, model-independent, estimates of the masses of the 
four components of this system. 
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